Novel exponents control the quasi-deterministic limit of the extinction transition 
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The quasi-deterministic limit of the generic extinction transition is considered within the frame- 
work of standard epidemiological models. The susceptible-infected-susceptible (SIS) model is known 
to exhibit a transition from extinction to spreading, as the infectivity is increased, described by the 
directed percolation equivalence class. We find that the distance from the transition point, and the 
prefactor controlling the divergence of the (perpendicular) correlation length, both scale with the 
local population size, TV, with two novel universal exponents. Different exponents characterize the 
large TV behavior of the susceptible-infected-recovered (SIR) model, which belongs to the dynamic 
percolation class. Extensive numerical studies in a range of systems lead to the conjecture that 
these characteristics are generic and may be used in order to classify the high density limit of any 
stochastic process on the edge of extinction. 
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INTRODUCTION 

Perhaps the most useful approximation in science is 
the description of microscopic, stochastic processes via 
"mean field" deterministic equations. Except for the 
most fundamental parts of particle physics, almost any 
other study involves that procedure: classical mechanics 
is used to describe quantum systems, classical electro- 
magnetism neglects the fluctuations of the photon den- 
sity and continuum mechanics averages out the micro- 
scopic stochasticity involved in the motion of an individ- 
ual molecule. Chemical reaction kinetics is generally de- 
scribed by rate equations; population dynamics and other 
ecological processes are depicted using the concept of lo- 
gistic growth, or by Lotka-Volterra equations. The dy- 
namics of the individual, microscopic, constituents is al- 
ways stochastic, and is subject to fluctuations; these fluc- 
tuations are smeared out when the system is described 
by deterministic rate equations. 

Basically, the underlying assumption beyond all these 
approximations is that the microscopic fluctuations are 
averaged out in the "large density" (many atoms, ani- 
mals, quanta) limit. A generic analysis of the determin- 
istic limit, like a 1/TV expansion where TV is the number 
of elements, is still lacking in many fields. Problems like 
the quantum classical correspondence in chaotic systems, 
or the decay of a quasistationary state to an absorbing 
state are still a subject of intensive studies. The situation 
becomes even more complicated when spatially extended 
systems, made of diffusively coupled patches, are con- 
sidered. What exactly determine the "large N" limit in 
that system? Should the number of microscopic entities 
be large on a single patch, or within a correlation length? 
All these questions are still open, as no systematic per- 
turbation theory in 1/TV exists so far. 

In this paper we consider a generic process in popu- 
lation dynamics: the extinction transition with an ab- 
sorbing state. The framework used is two well-known 
models for epidemics, namely, the SIS and the SIR mod- 



els [H, 0]. For a well mixed population of size TV (say, 
on a single patch) the stochastic process starts by the 
introduction of a single infected ( "I" ) individual into the 
system. This individual may infect any other, suscepti- 
ble ("S"), individual with rate a/TV, where after the in- 
fection the susceptible becomes infected and may infect 
other susceptible members of the community; the only 
other process is a "recovery" of an infected person; this 
happens with rate (3. In the SIR model, the recovered 
individuals are then immune against the disease, while 
the SIS model describes the case where the recovered be- 
came susceptible again. The elementary processes, thus, 
are 
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Denoting by / the fraction of infected individuals, and by 
S the fraction of susceptibles, the mean field equations 
in the well mixed limit are 
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Clearly, the SIR process is self-limiting, as S decreases 
with time, while the SIS process may support an endemic 
state with / = 1 — a/ (3. Another piece of information gar- 
nered from the rate equations is the existence of a tran- 
sition when the infectivity parameter Rq = a/ (3 crosses 
the value R c = 1. Below Rq — R c , both SIS and SIR pro- 
cesses immediately decay; above that value, an outbreak 
is possible. 

The mean field equations ([2|) are an approximation for 
the real stochastic process ([1]). Some characteristics of 
the epidemic, though, may be calculated exactly for the 
real stochastic process in the well mixed limit. In partic- 
ular, 

• Below R c both processes are the same in the large 
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TV limit, since the change of a susceptible becoming rein- 
fected in the SIS model is vanishingly small. 

• For large N, above R c there are two peaks in the 
distribution of epidemic sizes. The first peak is at size 
1 (i.e. the initial infected individual, only). The second 
peak for the SIR model is at the size predicted by the de- 
terministic equations For the SIS model, the second 
peak is at an exponentially large value, corresponding to 
the exponentially long lifetime of the metastable endemic 
state[i]. 

• At R c , the average size of the epidemic scales for 
large N as N 1 / 2 for the SIS model @,[| and as N 1 / 3 for 
the SIR model 0,0, S|. 

Let us consider now a one dimensional array of L 
patches, with N susceptible individuals on each patch, 
where a single infected individual is introduced at a sin- 
gle patch. 

Above R c , and for N — > oo, the deterministic equations 
(|2"|), with an appropriate contact term playing the role 
of "diffusion" , give rise to a constant velocity solution. 
The mean field dynamics is akin to that of the FKPP 
equation 0, [1(3] , and yield a velocity proportional to the 
square root of the distance from R c . The solution for the 
SIS case is that of a front separating the metastable un- 
infected state ahead of the front from the endemic state 
behind. In the SIR case, the solution is a pulse, leaving 
behind a state with a diminished susceptible population 
such that the effective infectivity parameter, RqS/N, is 
below R c . 

What happens in the stochastic models then? Let us 
discuss the SIS case first. For N — 1 the SIS model 



is equivalent to the contact process |ll|, where there is 
a transition to a propagating state at a finite value of 
Rq > 1. This transition is known to be in the directed 
peculation (DP) equivalence class [13]. Above the tran- 
sition, there is a finite chance of generating a wave of 
infection whose lifetime is infinite. The average velocity 
of the wave is less than that predicted by the determin- 
istic equation. This has been analyzed by Brunet and 
Derrida [HI], who show that as N — > oo, the average 
velocity approaches the deterministic value, albeit with 
anomalously large 0(ln~ 2 N) corrections. Below the di- 
rected percolation transition, the typical spatial extent 
of the epidemic, denoted £x> is finite, and diverges as the 
transition point is approached. Following the Janssen- 
Grassberger conjecture [3], it is widely believed that any 
generic extinction transition, and, in particular, the tran- 
sition for the SIS model for any N, falls in the equivalence 
class of directed percolation. 

In order to test these predictions in the large density 
limit we have simulated the SIS process on a one dimen- 
sional array of sites. On each site there are N individ- 
uals, and this number is fixed throughout the process. 
The epidemics is ignited by the introduction of a single 
infected individual on a single site. The chance of an 
/ to infect any susceptible individual on the same site 



is (1 — x)a/A r , and its chance to infect an S on one of 
the two neighboring sites is \ a /N. Thus, \ is the inter- 
site infectivity rate, while the population within a site is 
considered as well mixed. The chance of recovery is inde- 
pendent of the spatial structure and the recovery rate is 
(3 for any individual. We have used an exact, agent-based 
simulation in order to study the approach to transition 
from below, using the divergence of the perpendicular 
correlation length £j_ as a marker of the transition. 

What is the effect of increasing N on the transition? 
Fixing (3 and x one should expect a dependence of Rdp, 
the value of Rq at the DP transition, on N. If N = 1, for 
example, the SIS becomes a simple contact process on a 
line, and the DP transition happens when the infection 
rate is 3.297 times higher than the recovery rate; accord- 
ingly, R DP (N = 1) = 3.297/x- On the other hand, as 
N — ► oo, Rdp should approach the value R c = 1 and 
becomes x independent, as a single site may support the 
endemic state. Rdp should interpolate between these 
two limits as a function of N. 

Another effect has to do with the divergence of £j_ at 
the transition. According to the stochastic theory, right 
below the transition £j_ is finite but large. On the other 
hand, the classical equations predict that £j_ = below 
the transition. One expects both statements to be true 
at the large N limit. The only way out of that paradox 
is to understand that the classical description must fail 
close enough to the transition, but the region in which 
it fails must shrink to zero in the N — ► oo limit. Say it 
another way, at a fixed distance from the transition, and 
for large enough N, the perpendicular correlation length 
should approach zero. 

Our simulations support all parts of this picture. First, 
the Jansen-Grassberger conjecture works and in all cases 
studied we find that £j_ diverges as the transition point is 
approached from below with the expected power — v± w 
— 1.097. Second, we find that £j_ obeys a universal scaling 
with N, independent of the details of the model. Specif- 
ically, near the transition 



e± = A(N)(R DP (N) - Ro)- Vx 
where for large N 



(3) 



R DP (N)-R C ~N- K , k«0.66 (4) 
where for large N 



A(N) ~ N- 



0.41. 



(5) 



The data supporting these findings are presented in Fig. 
[T] for the specific case of an SIS model with fixed recov- 
ery time and x = 0.2. The applicability of these re- 
sults, though, is much wider; the same exponents turn 
out to describe the large N limit of the SIS model with 
different x an d with exponentially distributed recovery 
time, as well as other stochastic models like branching- 
annihilation-death (A — ► 0, A + A — > 0, A — > 2 A) and 
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FIG. 1: Directed percolation transition point, Rdp, of the ID 
SIS model as a function of TV, together with the critical per- 
pendicular correlation length amplitude A(N). Also shown 
are power-law fits to the data. £ = 0.2 



so on. We thus conjecture that these exponents describe 
the large TV behavior of any one-dimensional extinction 
transition that belongs to the DP equivalence class. 

The data collapse presented in Figure [2] reveals the 
existence of much stronger regularity at large TV's. Eq. 
(J31) implies that j\T re ~ T / I '- L £7 1 ' !/ ~ L is a linear function of 
(i? — R C )N K close to the transition. In fact, this result 
generalizes to a whole scaling regime where the scaled 
correlation length is a function of the scaled distance from 
the classical transition: 



N K - 



F((R - R C )N K ) 



(6) 



Again, we have verified that this scaling behavior is in- 
dependent of the (nonzero) strength of the intersite cou- 
pling. 

We now turn to examine the behavior of the stochastic 
SIR model at large TV. Here, the model is supposed to 
be in the dynamic percolation universality class 15, 
This class does not have a transition in one spatial dimen- 
sion, and the propagating pulse always dies out in an in- 
finite system. Still, the correlation length exhibits a scal- 
ing behavior with TV very similar to that of the SIS model. 
We find that the scaled correlation length, £^TV -0 - 209 is 
a function of the scaled infectivity (R — R C )N ' . The 
data collapse is presented in Fig. [3J Here the correlation 
length is finite for all infectivities, indicating that we are 
indeed below the percolation transition. The correlation 
length does however grow very rapidly with R a . 

In the simulations presented above, we have used as 
our measure of the correlation length the inverse of the 
exponential falloff rate of the distribution function for in- 
fection location. It is remarkable that in both models, for 
all R studied (below the extinction transition), this dis- 
tribution function appears to be given exactly (at least to 
within our statistical errors, over some seven decades) by 
the functional form P(x) = C{R 0l N)e~\ x - x °\/^ R °' N \ 
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FIG. 2: Collapse of the ID SIS model scaled correlation 
length, £J 1 ^" ± TV 0,284 versus the scaled infectivity (R a — 
l)iV - 662 for TV = 200, 400, 800, 1600 and 3200. The con- 
tact parameter is again \ = 0.2, and the fixed recovery time 
variant was used. 
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FIG. 3: Collapse of the SIR scaled correlation length, 
£ x TV-°- 209 versus the scaled infectivity {R a - 1)TV ' 454 for 
TV = 200, 400, and 800. Here too \ = 0-2. 



with x is the location of the original infected person. 
The average total number of infections is given then 
roughly by Af = 2C(i? , N)£ X (R , TV) (assuming £i » 1, 
as it is for TV large). We can then deduce the scaling prop- 
erties of C(R = 1, TV) with TV. The number of infections 
at the central site, C(1,TV) are the result of the on-site 
infections nucleated by the original infected person and 
by those nucleated by the neighboring sites. The effec- 
tive on-site infection rate is R (l — x), and since this is 
below R c = 1, each initial infection develops essentially 
independently, giving rise to 1/(1 — R (l — x)) total in- 
fections. The number of primary infections induced by 
the neighbors is xCexp(— l/£j_) ~ xC(l ~ V£-0- Thus, 
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we have, for R Q = 1 



C 



1 - 
= C4 

This implies that 



(1 
1 

X 



X) 

c_ 



c = 



1 + xc 1 



(7) 
(8) 

(9) 



so C(l, -/V) scales exactly the same with JV as £_i_(-ff = 1). 
Thus, the total number of infections scales as ti\(R = 
1), in both the SIS and SIR models. For the SIS 
model, this gives the total number of infections scaling as 
N 2(ku ± -t) _ ^0.623 ^ wnereas i n tne SIR model, we have 

the total number of infections scaling as approximately 
iV 418 . Notice that in both cases, the total number of 
infections at the critical infectivity scales faster than the 
dimensional results, 1/2 and 1/3, respectively. We 
have tested these predictions in our simulations, (data 
not shown) and found them to be very well satisfied. 
Also, preliminary work indicates that these exponents 
are higher still in two dimensions, and saturate at 1 as 
the dimension goes to infinity [l7j ]. 

In the SIS model, as in the related contact process, 
the total number of infections diverges at the directed 
percolation transition, with the exponent 7 = v\\ + v±_ — 
2/3 = 2.278. However, the total number of infections is 
equal to C£j_- Thus, C should diverge at the transition 
as 7 — u± = 1.181. Combining this with our above result 
for the scaling of C with N at criticality, we get that 



(CW 



-0.311-1-1/1.181 



g((R-i)N 



0.668-1 



(10) 



where Q vanishes linearly at the transition. It is remark- 
able that even though the directed percolation exponent 
associated with C differs from that of £j_ , the N scaling 
is the same. 

Our simulations suggest that the quasi-deterministic 
region is controlled by the new critical exponents k and 
t. At the transition, the deterministic limit does not ex- 
ist and the stochastic dynamics of the microscopic con- 
stituents determines the system behavior for any time 
scale. Off transition, on the other hand, the effect of 
stochasticity vanishes for large enough TV. The value of 
k and r is parameter independent, as long as the di- 
mensionality of the system and the type of stochastic 
transition (DP vs. dynamic percolation) are kept fixed. 

It is interesting to point out an exception to our anal- 
ysis, namely, the branching-annihilation process with no 
bare death term (A —* 2 A, A + A — > 0) for Brownian par- 
ticles. Here at low densities there is a DP type transition 
at finite birth rate, but as N — ► 00 (i.e., A — » 0) the tran- 
sition approaches a = 0. The deterministic limit of the 
transition in that case does not belong to the directed 



percolation equivalence class; in fact, it is known that 
for fj = the density decays like i -1 / 2 in one spatial di- 
mension. Indeed, Cardy and Tauber [l8| calculated that 
the distance to the DP transition should scale as 7V~ 2 , 
(since close to the classical transition the density is in- 
versly proportional to 1 /N) . Addition of a spontaneous 
death process A A to the model shifts the determinis- 
tic transition to a — fi and provide us with a model that 
admits a DP transition all the way to N = 00. In fact, 
our preliminary numerics show that in the Cardy- Tauber 
case the scaling of the critical a appears to be consistent 
with their prediction. Adding fi changes the picture dra- 
matically and o~ c ~ [i + N~ K , where k is identical to that 
measured for the SIS model. It is interesting to note that, 
if the Cardy- Tauber perturbative technique is modified 
to include the spontaneous death process, at large N the 
distance from the DP transition is predicted to scale as 
A^ 1 instead of A~ 2 . 

Many aspects of that problem are still open. In par- 
ticular a rigorous classification scheme for the quasi- 
deterministic behavior is still missing. The extinction 
transition in higher dimensions (where dynamic percola- 
tion admits a nontrivial critical point) and the behavior 
of the system above the transition also need to be in- 
vestigated. We hope to address these subjects in future 
work. 
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